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INTRODUCTION 

Monitoring of snow cover extent and snow water equivalent (SWT,) in boreal forests is 
important for determining the amount of potential runoff and beginning date of snowmelt. The 
great expanse of the boreal forest necessitates the use of satellite measurements to monitor snow 
cover (Derksen, 2008). Snow cover in the boreal forest can be mapped with either the Moderate 
Resolution Imaging Spectroradiometer (MODIS) or the Advanced Microwave Scanning 
Radiometer for EOS (AMSR-E) microwave instrument. The extent of snow cover is estimated 
from the MODIS data and SWE is estimated from the AMSR-E. Environmental limitations affect 
both sensors in different ways to limit their ability to detect snow in some situations. Forest 
density, snow wetness, and snow depth are factors that limit the effectiveness of both sensors for 
snow detection. Cloud cover is a significant hindrance to monitoring snow cover extent using 
MODIS but is not a hindrance to the use of the AMSR-E, These limitations could be mitigated by 
combining MODIS and AMSR-E data to allow for improved interpretation of snow cover extent 
and SWE on a daily basis and provide temporal continuity of snow’ mapping across the boreal 
forest regions in Canada. The purpose of this study is to investigate if temporal monitoring of 
snow cover using a combination of MODIS and AMSR-E data could yield a better interpretation 
of changing snow cover conditions. 

The MODIS snow mapping algorithm is based on snow detection using the Normalized 
Difference Snow Index (NDSI) and the Normalized Difference Vegetation Index (NDVI) to 
enhance snow' detection in dense vegetation (Riggs et al„ 2006). (Other spectral threshold tests are 
also used to map snow using MODIS.) Snow cover under a forest canopy may have an effect on 
the NDVI thus we use the NDVI in snow detection (Klein et al., 1998). A MODIS snow fraction 
product is also generated but not used in this study. In this study the NDSI and NDVI components 
of the snow mapping algorithm were calculated and analyzed to determine how they changed 
through the seasons. 

A blended snow product, the Air Force Weather Agency and NASA (ANSA) snow algorithm 
and product has recently been developed (Foster ct al., in press). The ANSA algorithm blends the 
MODIS snow cover and AMSR-E SWE products into a single snow product that has been shown 
to improve the performance of snow' cover mapping (Hall el al., 2007). In this study components 
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of the ANSA snow algorithm are used along with additional MODIS data to monitor daily 
changes in snow cover over the period of 1 February to 30 June 2008. 


STUDY AREA & DATA ANALYSIS 

A region of the boreal forest across northern Manitoba and Saskatchewan, Canada, 
fw\vw.mb,ec.gc.ca/'n ature/ecosvstems/wbci-icpo/di00s03.en.html) was chosen for this study. The 
daily MODIS data products are mapped in Sinusoidal projection and stored in tiles that cover an 
approximately 10° x 10° area of the earth surface. A region of approximately 970 x 600 km 
(Figures 1-3) was chosen for study. This region is covered by the MODIS daily products stored in 
the tile hl2v03 (horizontal and vertical tile indicies). MODIS snow cover and MODIS daily 
surface reflectance data products, and the AMSR-E daily SWE product were analyzed. The 
AMSR-E SWE daily product is global coverage in EASE-Grid projection; the study area was 
located by geographic coordinates in the products. The MODIS snow-covered area daily maps, 
MODTS surface reflectance products and AMSR-E SWE data product are available to monitor 
snow cover over the season but are used separately. The time period monitored for the study was 
from 1 February to 30 June 2008. The objective of the study was to determine if there are 
relationships among the snow products or spectral ratios of NDSI and NDVI that could be utilized 
to monitor snow cover, detect the beginning of snowmelt, and to interpret probable transient wet 
snow events that impact the quality of the SWE estimates. 



Figure 1 . View of the boreal study area shown as MODIS bands 1 , 4, 6 RGB image to highlight the 
presence of snow. Snow appears yellow in this band combination, pixels dominated by snow are brightest 
yellow; snow mixed with other surface features appears in decreasing shades of yellow. Clouds appear in 
white and light gray over some surface features. The study area is outlined (black polygon). Date of image is 
10 April 2008; the daily surface reflectance product MOD09GA, tile hl2v03, 500 m resolution, mapped in 
geographic projection. The star shows the approximate location of Thompson Zoo, Manitoba weather station. 

MODIS data products used are the daily snow cover (MOD10A1), and the daily surface 
reflectance (MOD09GA) both at 500 m resolution. Daily data were extracted for the boreal forest 
region and processed for each day of the study. The NDSI = (MODIS band 6 - band 4 / band 6 + 
band 4); NDVI = (MODIS band 2 - band 1 / band 2 + band 1). (MODIS central wavelength 



locations; band 1, 0.64 p; band 2, 0.85pm; band 4, 0.55 pm; band 6, 1.6 pm.) MODIS 
reflectance data from the MOD09GA daily surface reflectance product, 500 m resolution were 
used to calculate these indices. For each day, MODIS observations that were not cloud obscured 
were used for analysis. The cloud mask from MOD10A1, carried through from the MODIS cloud 
mask product in processing, was applied as the cloud mask. The number of observations used 
varied depending on the extent of cloud cover each day. Daily SWE data were extracted for the 
study region from the AMSR-E SWE product at 25 km resolution. 



Figure 2. MODIS snow map, 500 m resolution, of the boreal study area on 1 0 April 2008, tile hi 2v03 with 
the study area outlined (black polygon). Snow is white, snow and ice covered lakes are blue-green, clouds 
are pink. The image shown is mapped in geographic projection. The star shows the approximate location of 

Thompson Zoo, Manitoba weather station. 


DISCUSSION & RESULT 

Snow cover extent and SWE appear to remain relatively constant through the winter except for 
an event 8-15 April. After that snowmelt occurs over the period of about 20 April through early 
May as indicted by declining snow-covered area and SWE (Figure 3). Snow-covered area is 
100% through the winter though there is day to day variation in the NDSI and NDV1. The SWE 
amount is relatively constant through the winter though there are day to day variations. The day to 
day variations in snow-covered area and SWE are interpreted as normal and attributable to day to 
day changes in environment and observing conditions. Drastic changes over a few days are of 
interest as that signals a significant change in snow pack conditions or the start of snowmelt. 




Figure 3. Average daily values for MODIS snow-covered area, NDSI, NDVI and AMSR-E SWE in the 
boreal forest study region. Averages for the MODIS data were derived using only cloud free pixel 
observations of land surface. A possible wet snow event occurred between the dotted vertical lines. Figures 
1 , 2, 3B and 4 are from 10 April, noted as the dash-dot line on plot, a relatively cloudless day over the region. 

For the purpose of monitoring the snow cover detecting the start date of snowmelt is important 
as the hydrologically significant SWE is determined at that time. The start of snowmelt appears to 
be detected by the MODIS snow-covered area, NDSI, and NDVT showing a rapid decreasing 
trend in snow-covered area beginning about 20 April (Fig. 3). Decreases in the snow-covered area 
and NDSI with corresponding increase in NDVI signal the beginning of snowmelt on about 20 
April (Fig. 3). Those indicators presage a rapid decreases in SWE beginning about 1 May (Fig. 3). 
The SWE also begins a decreasing trend a few days after the MODIS snow cover indicators began 
decreasing. However, a very rapid decrease in SWE followed by a rapid steep rise in SWE 
occurred several weeks before seasonal melt of the snow cover began. That anomaly in SWE 
trend was likely caused by a wet snow event. It is possible that the wet snow event could have 
been misinterpreted as the start of snowmelt if SWE monitoring alone was used. 

Rapid changes in the wetness of a snow pack can cause the snow pack to ‘disappear’ to the 
AMSR-E because wet snow has a negligible scattering signal (Kelly et al., 2003) but if the snow 
refreezes the scattering signal returns causing the snow to re-appear. A wet snow event appears 
to have occurred between 8-15 April, as seen by the dotted vertical lines on Fig. 3. Evidence for a 
possible wet snow event is the sudden rapid drop in derived SWE (increases in Tb) with a sudden 
sharp rise in SWE after a few days to SWE values in the range prior to the drop. SWE maps over 
the period 6-19 April (Fig. 4) chronicle the areal extent of the decline and rise in SWE. 
Meteorological data from the Thompson Zoo, Manitoba station listed in Fig. 5 indicate a warm 
period with continuous snow cover over that period. The MODIS snow-covered area, NDSI and 
NDVI remained relatively consistent over that period (Fig. 3) indicating continuous snow cover, 
though NDSI did spike during the period. Conditions observed on 10 April, early in the warm 
period, by MOOTS, AMSR-E and ANSA are shown in Figs. 1, 2, 4B and 6. Though some of the 
snow cover may have melted and depth decreased (Fig. 5), the snow cover did not almost 




disappear as the SWE maps show on 14 April. The MODIS snow cover map, NDSI, NDVI, and 
the ANSA map all show a continuous snow cover over the wet snow period. In the ANSA the 
snow cover is continuous and some amount of SWE (Fig. 6) is shown in the region. In the 
classification scheme used in the ANSA, SWE is not broken into ranges for combination with 
MODIS snow cover area. There was no reported snowfall (Fig. 5) that would subsequently cause 
the SWE to return to the range shown for about 6 April on 1 9 April. The rapid divergence of SWE 
compared to the MODIS snow cover indicators suggest that SWE estimates were affected by a 
transient wet snow event, which caused SWE estimates to fall to less than 50% of pre-event values 
then suddenly rebound. 



Figure 4. Sequence of AMSR-E SWE maps, 25 km resolution, of the boreal forest study area (black polygon) 
over the time of the probable wet snow event. (From left to right, A,B,C,D.) Dates shown arc A) April 6, B) 
April 10, C) April 14, D) 19 April 2008, mapped in geographic projection. The star shows approximate 
location of Thompson Zoo, Manitoba weather station. 




SWv -v. 



■ tf* ;&* 3e* m<*/ 




A • .?• r 

S:.w 






■ ■ 


mn 

ia-i“ 


P 


.... HAi 


• A.ri v” . 






i l r *iOa*Y OtwnNStion&at 

k | S3 1 Xew ta 



< PR?5^5S'22?t# 



A triages. 









A 


Canadian 

Dally Dsta Repo? 

for April 20QS 















. ^ ^ 

;. (CfXtVf 

Note* on Q£txJ3U8&i st. 








Cantjtct* *■ 



m&wfiovziX} 

















55“ 4S.0C0’ N 






warn 

Question# 



WWlfttiW 9'* 52.200 W 


205.70 m 




WMQtPl 



W.Br 



• unfc* * 







M f. i v 

1 Water Data 


J April 

~3!5o«3_ 

68 f 


asKJ^asli ■ 


1 


Daily D«?«s ter April OTO0 


1 


I PrcaaJy® Qisdeaur* 

tt 

Site 

Mm Tmm Xwm 


goat im 

Idtffi 

®LBt &MM 1 



4 

!>%& 

«*u Safes gww 

t'9 HSK 


,;i j-';: t 


y 

°-n. *e •'C 

SMk 

dtoS wn o» 

w 

am) urn 

IMtU 

' -v "M 



Id afi 

■-7C- . 

" »C . M 

d 

*£ 

cm JO’s 

km/h 





fel 

a; 







tut 

3.0 -1/0 -7.0 

2S.0 

0.0 0.0 

0.1 

0 1 

ft/ 




t)3?f 

5,0 -9.D -2-0 

20.0 

0.0 0.0 

0:9 

03 

66 





•13 5 -21,5 -17,5 

35.5 

0.0 D.O 

D.O 

DO 

fti- 


. ■' V' 

> 

04 1 

.••16,5 -30.0 -22.0 

40.3 

0,0 5,0 ■ 

M 

02 

65 




«&♦ 

-9.S -20.5 -15,0 

33.0 

0,0 0.0 

0.0 

0 0 

65 





Of. I 

0.0 -23.0 -U.S 

39. S 

o.o s.n 

0.0 

0.0 

6S . 


■ 

|| t"' 


a/s 

5,5 -23.5 -9.0 

27.0 

0.0 0,0 

a.o 

O.D 

65 




| 

mt 

6,0 -SO. 5 -2,3 

20.3 

045 0,0 

3,0 

0.0 

64 





mi 

7.S -7.0 03 

17.7 

0.0 0.0 

3,0 

t».u 

62 




i 

Hit 

ta.s -2,« 

14.0 

OX! 0.0 

s.o 

0,0 

se 





5 S t 

6.S -l.S 2.5 

15,5 

0,0 o.c 

0,0 

0.0 

S7 





12* 

10,0 -2.5 3.6 

l*,2 

0,0 8,0 

0.0 

0.0 

£4 



•Sill 


i*t 

1S.S 3.S 9.S 

8,5 

0.0 0,3 

a •• 

o.e 

47. 



IliilM 


is* 

I2.fi -2,0 S.3 

15.7 

D.O G.S 

0.0 

0.0 

40 





16* 

. 3.5 -3,0 1.5 

• 16,7 

3.0 3,0 

0,0 

o.c 

30 



■fp 1 


f7t 

7.G -7-S -0,3 

18.3 

D.O D.O 

0,0 

o.c 

30 





t&i 

14,0 0,5 7,3 

10.7 

0,0 . o.c . 

043 

0.0 

38 



" , ^ /j, J 


5:9 f 

4,6 -2.0 1,3 

16,7 

0.0 G.G 

00 

0.0 

26 





nm 

0.0 -7.S -3s» 

21.8 

6-0 0,0 

0,0 

0,0 

sb 





21* 

2.S -4.0 -0.8 

18,8 

0,0 0:0 

1,7 

1.7 

22 





22* 

8.0 -4,5 -1.8 

13.2 

8.0 GO 

0:0 

0.0 

St 





3ST 

-l.S -55 -3.S 

S3,s 

0.0 0 0 

0.3 

0.3 

IS 





24 i 

2.6 -3*43 -6,0 

24,0 

0,0 0 0 

0,0 

00 

19 





25 1 

2.5 -10.11 -3.8 

21.8 

e,e- o.o 

on 

6.D 

16 





2&r 

3,3 -11.0 -3,8 

25,8 

o.o o,a 

D:S 

•o.e 

IS 



: I 


27 f 

3.0 -10.0 -3,6 

21.5 

0.0 o.o 

O.fJ 

0 0 

16 



" 'y ■ , J, j* 

i&ena 

..,WKS~. 









j 


j ^ AW Stag party j __ . j j 


Figure 5. Daily data report April 2008, Thompson Zoo, Manitoba. 
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(201) MODIS snow 1-100% and SWE land not processed 
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(0) Land 
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(1498) Fill 


Figure 6. ANSA snow map, 25 km resolution, of the boreal forest study area on 1 0 April 2008, shown 
mapped in Lambert Azimuthal projection. The star shows approximate location of the Thompson Zoo 

Manitoba, station. 


The suggestion emerging from this study is that monitoring and analysis of trends in snow- 
covered area and other indicators of snow cover observed concurrently with MODIS and AMSR-E 
could be utilized to interpret transient wet snow events that affect SWE estimates. In addition, 
trends of these snow cover indicators could be used to determine when seasonal snowmelt begins. 
An anticipated benefit would be improved estimates of the extent, volume and timing of 
snowmelt. 
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